Level spacing distribution of critical random matrix ensembles 
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We consider unitary invariant random matrix ensembles which obey spectral statistics different 
from the Wigner-Dyson, including unitary ensembles with slowly (~ log 2 x) growing potentials and 
the finite-temperature fermi gas model. If the deformation parameters in these matrix ensembles are 
small, the asymptotically translational-invariant region in the spectral bulk is universally governed 
by a one-parameter generalization of the sine kernel. We provide an analytic expression for the 
distribution of the eigenvalue spacings of this universal asymptotic kernel, which is a hybrid of 
the Wigner-Dyson and the Poisson distributions, by determining the Fredholm determinant of the 
universal kernel in terms of a Painleve VI transcendental function. 

PACS numbers: 05.45. +b, 05.40. +j, 71.30.+h 



I 

C/3 



I 

C 

O 
o 



> 

00 
(N 
l> 
O 
oo 

C3 



T3 

C 

O 
O 



X 



A main goal of quantum chaos study is to describe 
quantitatively statistical behaviors of spectra of classi- 
cally non-integrable systems, such as complex nuclei [Q, 
billiards Q], QCD ||, and disordered systems [Q. A char- 
acteristic observable in such studies, used analytically or 
numerically to measure the deviation from integrability, 
is the probability E(s) of having no energy levels in an 
interval of width s, or the distribution of spacings be- 
tween adjacent levels P(s) — E"(s). These observables 
capture the behavior of local correlations of large number 
of energy levels, as the former consists of an infinite sum 
of integrals of regulated^ spectral correlators, 



~ f-l)p r s / 2 

s ( s )=ES l / dx ^ 



■dx p (p(xx) ■ ■ ■ p(x p )) . (I) 



A technical virtue of invariant random matrix models of 
quantum chaotic systems HQ is that any p-point spec- 
tral correlation function of the former can be composed 
from the connected two-point function as 



{p(xi)---p(x p )) I 



det K{xi,Xj). 



(2) 



This can most easily be proved by the orthogonal polyno- 
mial method || . It allows a compact expression for the 
level-free probability E(s) as a Fredholm determinant, 



E(s) = det(l - K) 



(3) 



over the intervall— s/ 2, s/2]. 

Jimbo et. al. |]7j have made a remarkable observation 
that the logarithmic derivative of the Fredholm determi- 
nant, R(s) = — (log -E(s))', of the sine kernel 



K(x,y) = 



sin7r(a; — y) 
ir(x - y) 



(4) 



describing the bulk correlation of the Gaussian unitary 
ensemble, satisfies the cr-form of a Painleve V equation, 

(#(«) + -#'(«)) 2 +(ir S R'{ S )) 2 = R'{ S ){R{s) + S R'{ S )) 2 {5) 



Their method is subsequently generalized by Tracy and 
Widom S to kernels of the form 



K{x,y) = 



(j>{x)ip{y) - ip(x)<j>(y) 



(0) 



whose component functions satisfy a set of first order lin- 
ear differential equations with meromorphic coefficients. 
This class of kernels includes AiryJ |k Bessel and 
their multicritical generalizations Jll Jl2| , describing cor- 
relations at the edges of the spectral bands. 

However, such invariant random matrix ensembles are 
sometimes crude idealization of physical systems in con- 
cern, based only on the symmetries of the systems. It is 
a priori unclear that random matrix ensembles can still 
provide quantitative description of realistic physical sys- 
tems where localization of the states can occur If 
there exists such a random matrix ensemble, it must be a 
non-trivial deformation of the classical invariant random 
matrix ensembles, so as to violate the wide universal- 
ity that the Gaussian ensembles possess |Q. One such 
example is a random banded matrix ensemble fill , mod- 
eling quasi ID materials. Another example is a random 
hamiltonian consisting of a Gaussian random matrix H 
and a diagonal real random matrix V, H( a ) — H+aV. In 
such cases spectral correlation functions generally do not 
allow the determinant form (|^), and E(s) is usually evalu- 
ated only perturbatively in s, by computing each p-point 
correlation function 119]. However this is not sufficient 
to determine P(s) = E/'(s) for large s, since it typically 
takes the form P(s) ~ s a exp(— cst.s 6 ) (generalizations 
of the Wigner surmise by Brody [flTj and by Berry and 
Robnik |18| ), whose exponential damping is invisible in 
the small-s expansion. Alternatively, by treating evo- 
lution in a as a diffusion process the joint probability 
distribution can be derived jl9|], but the level spacing 
distribution is yet to be obtained. The aim of this Letter 
is to derive P(s) from random matrix ensembles which 
describes a deformation of the Wigner-Dyson statistics 
while preserving unitary invariance at the level of parti- 
tion function. 
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Muttalib et. al. pQ] have introduced the g-analogue of 
the GUE with a potential (0 < q < 1) 

oo 

V(X) = ^log(l + 2 9 "cosh(2arcsinhA) +q 2n ). (7) 

Tl=l 

They have shown that, after unfolding the spectrum 



Ahi = J p(X)dX, 



this ensemble is described by a kernel 



K(x,y) = 



f(x + y) 6>i(7r(x-y),e " ) 
y/f(2x)f(2y) sinha(x-y) 

#4(7TX, e~~) 



(8) 



(9) 



cosh air 



1, 1 
a = - log-. 

2 g 



For e * 2 l a <C 1, there exists an asymptotically transla- 
tional invariant region where ([)]) is approximated by 



K(x,y) 



asin7r(a; — y) 
7rsinha(x — y) 



(10) 



It is clear from this form of the kernel that a set of eigen- 
values with \xi — Xj \ <C 1/a obeys the Wigner statistics, 
and that with |xj — Xj\ ^> 1/a obeys the Poisson statis- 
tics, i.e. is uncorrelated. Canali and Kravtsov [^lj argued 
that the U(N) symmetry is spontaneous broken in this 
ensemble, which induces a preferred basis in the matri- 
ces and deforms the statistics. The universality of this 
asymptotic kernel (^) observed for the g-Laguerre uni- 
tary ensemble [ p2| and subsequently proved for ensembles 
with potentials which grow very slowly as 



y(A)~^(log|A|) 2 (A -oo). 
2a 



(11) 



|21[| . This universality can be considered as an extension 
of Brezin and Zee's universality |l4| of the sine kernel 
for polynomially increasing potentials, who proved it by 
deriving the asymptotic form of the wave functions 



V>tv(A) ~ cos 7r / p(X)dX + 



if (A, A 



N-k \ 



M sin(7r(/ A p - J A p)) 
A - A' 



(12) 
(13) 



For polynomially increasing potentials, the spectral den- 
sity is bounded and is locally approximated by a con- 
stant. Therefore the unfolding is just a linear transfor- 
mation, leading universally to the sine kernel (|J). How- 
ever, for the potential ([Tl|), the spectral density behaves 
as p(A) ~ l/(2aA) implying an unusual unfolding map 
A — x = l/(2a)logA, while the formula ( |l3| ) stays valid 
p3[ . Then the kernel ( |l3| ) universally reduces to ( |l(i| ) 
after this unfolding. 



Chen and Muttalib [£4| have interpreted a particular 
unitary ensemble with V(x) ~ (log a;) 2 as a fermionic sys- 
tem at finite temperature. This link is made more con- 
crete by Moshe, Shapiro, and Neuberger |^5| who have 
introduced a random matrix ensemble 



Z= I d N He~ UH 



dUe -btr[U,H]lU,H]\ (14) 

U(N) 

as yet another unitary invariant deformed ensem- 
ble. For a given unitary matrix [/, the interaction 
btr[U, H][U, Hy tries to align random hermitian matri- 
ces H so that [U, H] = 0. The integration over U then 
amounts to recovering the U(N) invariance of the model 
which the GUE has enjoyed. The basis preference is still 
realized through the spontaneous breaking of the U(N) 
invariance. After integrating over U, the above model is 
identical to a system of ID free nonrelativistic fermions 
in a harmonic potential with curvature m = y/l + 4t> 
and at finite temperature (3 = arccosh(l + 1/26), extend- 
ing the well known fact that the GUE is equivalent to 
free harmonic fermions at zero temperature. Note that 
p-point correlation functions of the local spectral densi- 
ties p(X) = J2n=o l^n(A)| 2 /(e /3(c '*" M) + 1) are thus still 
expressed as determinants of the kernel 



K(X,X') = Y1 



^n(A)C(AQ 



(15) 



Accordingly its level-free probability can be expressed 
in terms of the Fredholm determinant of (|l5|). Using 
the asymptotics of the one-particle wave function ipn{x) 
given by the Hermite polynomials, these authors have 
also obtained the local form of the kernel (|l(i| ) with 
a = ir 2 /(2N(3). There the limit N — > oo is taken while 
keeping the microscopic unfolded variable x fixed. The 
formula ( |l0| ) is valid as long as N(3 is not too high to 
invalidate the grand canonical picture. The models jll| ) 
and ( |l4| ) are subsequently unified as an ensemble with 
multifractal eigenvectors p6[ , and the parameter a is 
identified as a measure of the multifractality. 

Surprisingly, this universality within random matrix 
theories, if extended to orthogonal ensembles (27), is ob- 
served to encompass the 3D Anderson model, i.e. a 
particle hopping on the lattice with random cite ener- 
gies. Canali fim has compared his Monte Carlo results 
of P(s) for the orthogonal ensembles with (l/2a) log 2 a; 
potentials with that for the Anderson Hamiltonian at 
the metal/insulator transition measured precisely in ref. 
p9[ by exact diagonalization, and observed an excellent 
agreement by tuning the coefficient to a « 2.5. There a is 
interpreted as the inverse dimcnsionless conductance at 
the transition point. Motivated by this success, we derive 
an analytic form of its level spacing distribution of ( 10) in 
this short Letter. Although the a — oo limit of the model 
does not obey the Poissonian statistics |50) as is naively 
expected, the error involved in the asymptotic kernel ( |To| ) 

is exponentially small (of order 0(e~ 7r '")) in the above 
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parameter range. We complete earlier attempts which 
computed P(s) numerically f^of or asymptotically 

We notice that the kernel ( |10[ ) is equivalent to that of 
Dyson's circular unitary ensemble at finite- N M: 



K(x,y) 



sin T (a; - y) 



(16) 



N sin \[x — y) 
by the following analytic continuation: 

N — , x^—x (17) 
a i 

Tracy and Widom H h ave also proved that the diagonal 
resolvent kernel of ^1^) is determined by a second-order 
differential equation which is reduced to a Painleve VI 
equation 32 1. Below we will reproduce their method. 
The kernel ([L0|) is written as 



K(x,y) = 



4>{x)ip{y) - ip(x)<f>(y) 



^2ax 



e 2ay 

ip{x) 



(18) 



These component functions satisfy 

<j> = a<j) + 7r^> , iji — — Tt(j) + aip. (19) 

We use the bra-ket notation <j>(x) = (x\(j>) and so on jj3|. 
Due to our choice of the component functions to be real 
valued (unlike @, sect. V.D), we have (x\6\<f>) = (4>\d\x) 
and similarly for ip with any self-adjoint operator O and 
real x. Then ( |l8| ) is equivalent to 

[e 2aS ,i,] = |0)(^|-|^)(0|, (20) 

where x and K are the multiplication operator of the 
independent variable and the integral operator with the 
kernel K(x,y)9(y — ti)9(t 2 — y), respectively. Below we 
will not explicitly write the dependence on the end points 
of the underlying interval [ti, t 2 \. The resolvent kernel is 
defined as 



R{x,y) = 
It follows from ( p(i|) that 

k , 1 



K 



1 - K 



\y)- 



l< 2 



MM] -MM): 
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1-K 1-K 1-K 
that is, 

(e 2a. _ e 2ay )R{ ^ y) = Q {x)P{y) _ P ( X )Q(y) } (23) 

Q(x) = (x\(l - K)- 1 ^), P(x) = (x\(l - k)- 1 ^). 
At a coincident point x — y we have 



2ae 2ax R{x,x) = Q'{x)P{x) - P'{x)Q{x). 



Now, by using the identity 
dk 

we obtain 



= (-l) i K\t i )(t i \, (i= 1,2) 



(24) 



(25) 



dQjx) 

dti 
dP(x) 



= (-iyR(x,t t )Q(ti), 

= {-iyR{x,u)p(u). 



(26a) 
(26b) 



On the other hand, by using the identity (D is the deriva- 
tion operator) 

[D,K]=K(\t 1 )(t 1 \-\h)(t 2 \), (27) 

which follows from the translational invariance of the ker- 
nel (d x + d y )K(x — y) = 0, we also have 

^. = {x \D(i-k)-^) 

= (x\(i - k)- 1 ^') + (x\(i - k)- l \D,k\{i - ky 1 ^) 

= aQ{x) + ttP{x) + R(x, ti)Q(ti) - R(x, t 2 )Q{t 2 ), (28a) 
dP{x) 
dx 

-%Q(x) + aP(x) + R(x, ii)P(ti) - R(x, t 2 )P(t 2 ). (28b) 

Now we set t\ = —t, t 2 = t, x,y = —t or i, and introduce 
notations q = Q(—t), q = Q{t),p = P(—t),p = P(t), and 
R = R(-t,t) = R(t,-t),R = R(t,t) = R(-t,-t). The 
last two equalities follow from the evenness of the kernel. 
Then eqs.@ and (§|) reads, after using @, 

(29a) 
(29b) 

(29c) 



pq — qp = 2R sinh 2at, 
2 ~ 

q 2 = - (R 2 sinh 2at + Ra e" 



2ul 



p~ + q 



IT 

2 



), 



(R 2 sinh2crf + Rae 2at ) 



The total i-derivatives of eqs. (|2^) lead to ( 
VP 

R = 2R 2 , R = ARR 



d/dt) 



qq = —(R sinh 2at) ' , 

7T 



(30) 



The left hand sides of eqs. (29) and 
tional constraint 



(31) 

satisfy an addi- 



(pp + qq) 2 + (pq - qp) 2 = (P 2 + q 2 ){p 2 + q 2 )- (32) 

By eliminating p,p, q, q, R and R from eqs.(p9|)-([32|), we 
finally obtain for R(s) (s = 2t,' = d/ds) 

2 



(21) 



(22) [a cosh as R'(s) + 



sinh< 



R"(s) +(tt sinh as R'(s) f 



R'(s) ((aR(s)) 2 + a sinh2as R(s)R'{s) + (sinh as R'(s)) 2 ). 

(33) 

This is our main result. It is equivalent to eq.(5.70) of 
ref. || after the analytic continuation ( p7| ) accompanied 
by a redefinition R(s) — > ^-R(s)} It clearly reduces to 
the Painleve V equation (0) for the GUE as a — > 0. In 
the Wigner-like region as< 1, we can expand hyperbolic 
functions into the Taylor series. Then (g) and (§ yield 

£?(«) = 1 - s + 0(s 4 ), 

E(s) = -(log^(s))' = 1 + s + 0{s 2 ), (34) 

By imposing this boundary condition, we obtain a per- 
turbative solution to eq.(33), 



R(s) =1 + s + s 2 



2 i 2 

7r +a 



5(7r 2 + g 2 ) \ 4 



P(s) = (i?(s) 2 -i?'(s))exp 



1 - 



(7r 2 +a 2 )(75-4^ 2 -6a 2 ) 
450 



7r 2 +a 2 2 (7r 2 +a 2 )(27T 2 + 3a 2 ) 4 (vr 2 + a 2 )(7r 2 + 2a 2 )(37r 2 + 5a 2 ) 6 (vr 2 + a 2 ) 2 (7r 2 + 4a 2 ) 7 

S — — S T" T~~ S — , ~ ~ ~ 5 



45 



945 



4050 



(35) 



,(36) 



which is in accord with the expansion (Q), (||). 

In Fig.l we show plots of the level spacing distributions 

P(s) for various a such that e~~ <C 1, obtained by nu- 
merically solving (33) under the boundary condition (34). 
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FIG. 1. Level spacing distributions P(s) of the kernel 
([[o|). The limiting distribution for a — > oo (ref. [30], eq.(85)) 
and the Poisson distribution are plotted for comparison. 

These distributions are indeed hybrids of the rescaled 
Wigner-Dyson distribution P(s) ~ s 13 (/3 = 2) for s^l/a 
and the Poisson distribution P(s) ~ e - cst s f or s^l/a. 
The extension of our result to the case of orthogonal 
ensembles (j3 = 1), which correspond to the Anderson 
model, will be reported elsewhere. 
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This is not trivial because ref. |8| has usedp = p* and q = ponent functions, cf>(— x) = cf>(x)* and ip(— x) — i/j(x)*. 

q* which follow from the analytic properties of its com- 
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